Explicit equilibrium modeling of transcription-factor binding and gene regulation

We have developed a computational model that predicts the probability of transcription factor binding to any site in the genome. GOMER (generalizable occupancy model of expression regulation) calculates binding probabilities on the basis of position weight matrices, and incorporates the effects of cooperativity and competition by explicit calculation of coupled binding equilibria. GOMER can be used to test hypotheses regarding gene regulation that build upon this physically principled prediction of protein-DNA binding

Antifungal drug resistance evoked via RNAi-dependent epimutations

Microorganisms evolve via mechanisms spanning sexual/parasexual reproduction, mutators, aneuploidy, Hsp90, and even prions. Mechanisms that may seem detrimental can be repurposed to generate diversity. Here we show the human fungal pathogen Mucor circinelloides develops spontaneous resistance to the antifungal drug FK506 (tacrolimus) via two distinct mechanisms. One involves Mendelian mutations that confer stable drug resistance; the other occurs via an epigenetic RNA interference (RNAi)-mediated pathway resulting in unstable drug resistance. The peptidyl-prolyl isomerase FKBP12 interacts with FK506 forming a complex that inhibits the protein phosphatase calcineurin1. Calcineurin inhibition by FK506 blocks M. circinelloides transition to hyphae and enforces yeast growth2. Mutations in the fkbA gene encoding FKBP12 or the calcineurin cnbR or cnaA genes confer FK506 resistance (FK506R) and restore hyphal growth. In parallel, RNAi is spontaneously triggered to silence the FKBP12 fkbA gene, giving rise to drug-resistant epimutants. FK506R epimutants readily reverted to the drug-sensitive wild-type (WT) phenotype when grown without drug. The establishment of these epimutants is accompanied by generation of abundant fkbA small RNA (sRNA) and requires the RNAi pathway as well as other factors that constrain or reverse the epimutant state. Silencing involves generation of a double-stranded RNA (dsRNA) trigger intermediate from the fkbA mature mRNA to produce antisense fkbA RNA. This study uncovers a novel epigenetic RNAi-based epimutation mechanism controlling phenotypic plasticity, with possible implications for antimicrobial drug resistance and RNAi-regulatory mechanisms in fungi and other eukaryotes

Environmental and Genetic Determinants of Colony Morphology in Yeast

Nutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs

Data from: The genetic architecture of biofilm formation in a clinical isolate of Saccharomyces cerevisiae

Biofilms are microbial communities that form on surfaces. They are the primary form of microbial growth in nature and can have detrimental impacts on human health. Some strains of the budding yeast Saccharomyces cerevisiae form colony biofilms, and there is substantial variation in colony architecture between biofilm-forming strains. To identify the genetic basis of biofilm variation, we developed a novel version of quantitative trait locus mapping, which leverages cryptic variation in a clinical isolate of S. cerevisiae. We mapped 13 loci linked to heterogeneity in biofilm architecture and identified the gene most closely associated with each locus. Of these candidate genes, six are members of the cyclic AMP-protein kinase A pathway, an evolutionarily conserved cell signaling network. Principal among these is CYR1, which encodes the enzyme that catalyzes production of cAMP. Through a combination of gene expression measurements, cell signaling assays, and gene overexpression, we determined the functional effects of allelic variation at CYR1. We found that increased pathway activity resulting from protein coding and expression variation of CYR1 enhances the formation of colony biofilms. Four other candidate genes encode kinases and transcription factors that are targets of this pathway. The protein products of several of these genes together regulate expression of the sixth candidate, FLO11, which encodes a cell adhesion protein. Our results indicate that epistatic interactions between alleles with both positive and negative effects on cyclic AMP-protein kinase A signaling underlie much of the architectural variation we observe in colony biofilms. They are also among the first to demonstrate genetic variation acting at multiple levels of an integrated signaling and regulatory network. Based on these results, we propose a mechanistic model that relates genetic variation to gene network function and phenotypic outcomes

Individual Segregant Phenotype and Genotype Table

This file compiles phenotype and genotype data collected on individual segregants from the Bulk Segregant Analysis. Each row in the file contains data on a single segregant

Whole-genome comparison of Leu3 binding in vitro and in vivo reveals the importance of nucleosome occupancy in target site selection

Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein–DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein–DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide

Bulk Genotype Table

This file contains the allele counts (number of occurrences of an allele in a sequencing read) for each heterozygous site found in YJM311 from the YJM311 (parental), simple, intermediate, and complex se- quencing runs, and the G-statistic values at each of these sites for the pairwise comparisons between the simple, intermediate, and complex pools. Note that heterozygous sites where neither allele matched the reference sequence are very rare, so were excluded from these results. Each row in the file contains data on a single heterozygous site identified in YJM311. The coordinates are based on Saccharomyces Genome Database (SGD) release r62, (release date: February 18, 2009)